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ABSTRACT 

We measure hot and cold spots on the microwave background associated with supercluster and supervoid 
structures identified in the Sloan Digital Sky Survey Luminous Red Galaxy catalog. The structures give a 
compelling visual imprint, with a mean temperature deviation of 9.6 ± 2.2/iK, i.e. above 4cr. We interpret this 
as a detection of the late-time Integrated Sachs-Wolfe (ISW) effect, in which cosmic acceleration from dark 
energy causes gravitational potentials to decay, heating or cooling photons passing through density crests or 
troughs. In a flat universe, the linear ISW effect is a direct signal of dark energy. 

Subject headings: cosmic microwave background — cosmology: observations — large-scale structure of uni- 
verse — methods: statistical 



L INTRODUCTION 

The cosmic microwave background (CMB) is a snapshot 
of the early universe; however, the light we observe has 
been processed by large-scale structure at low redshift, in 
part through the late-time integrated Sachs-Wolfe (ISW) ef- 
fect (jSachs & Wolfe 1967 ). As photons travel through time- 
varying gravitational potentials, they are slightly heated or 
cooled. In a universe dominated by dark energy, the gravi- 
tational potential decays with time even in linear theory, heat- 
ing photons traveling through crests and cooling photons in 
troughs of large-scale matter density fluctuations. Hereafter, 
'ISW' refers to the full nonlinear late-time ISW effect, also 
known as the Rees-Sciama effect ( |Rees & Sciama|1968| l. 

The ISW effect from dark energy can be detected with 
the cross-correlation function between the projected galaxy 
density and microwave background temperature over the sky 
( [Crittenden & Turok 1996 ). Measurements from individual 
galaxy surveys detect the effec t with signal-to-noise no higher 
flian 3 ([Scranton et al. 2003 ; B oughn & Crittenden|2004 



shordi et aL p004; Padmanab han et al.||2005t |Raccan"e 



h et 



al. 2008). Recently, various groups have combined multiple 
datasets to arrive at a detection as high as 4.5(7, though er- 
ror estimation with correlated galaxy datasets complicates the 
physical interpretation (IHo et al. 2008 [ Giannantonio et al. 
|2008| l. Additionally, studies using wavelet analyses have sug- 
gested that the signal can be localized to particular regions on 
the sky that depend on both the CMB and the galaxy density 
(|McEwen et al. 2008 ). 

The ISW signal peaks at spherical multipole 20 at 
z = 0.5 in the galaxy-CMB cross-power spectrum l{l + l)Ci 
(e.g. 'Padmanabhan et al. 2005| l. This corresponds to struc- 
tures with angular radius ~ 4°, or ^ 100 /i^^Mpc. We 
call these large structures 'supervoids' and 'superclusters,' but 
they may be thought of as gentle hills and valleys in the lin- 
ear density field. In a ACDM universe, the ISW signal from 
these broad, linear over- or under-dense structures is expected 
to dominate over smaller-scale fluctuations in the density. 

In this study, we identified a sample of supervoids and 
superclusters in a galaxy survey that could potentially pro- 
duce measurable ISW signals. We analyze these structures 
by stacking cutouts of the CMB centered on their projected 
locations. Atrio-Barandela et al. (2008) recently used a sim- 
ilar method on much smaller spatial scales, stacking WMAP 



data behind clusters to detect the frequency-dependent ther- 
mal Sunyaev-Zeldovich ( [Sunyaev & Zeldovich|1972| effect. 
We show that our structures are, on average, associated with a 
significant temperature imprint on the CMB. This is arguably 
the first visually compelling detection of the ISW effect. In 
our conclusions, we discuss the application of this work to 
dark energy and analysis of secondary CMB anisotropics. 

2. DATA AND METHODOLOGY 

We used a sample of 1.1 million Luminous Red Galax- 
ies (LRGs) from the Sloan Di gital Sky Survey (SDSS) 
( Adeknan-McCarthy et al. 2008) l covering 7500 square de- 
grees about the North Galactic pole. They span a redshift 
range of 0.4 < z < 0.75, with a median of ^ 0.5, and 
inhabit a volume of about 5 h^^G^c" . LRGs are elliptical 
galaxies in massive galaxy clusters representing large dark- 
matter halos (Blake et al. 2008), and are thought to be physi- 
cal ly similar objects across their redshift range ( Eisenstein et| 
|aL][200T] [Wake et al.|[2008| ). This makes them excellent, al- 
beit sparse, tracers of the cosmic matter distribution on scales 
> 10 Mpc. Our sample was selected from photometric data 
based on the criteria used in the Mega-Z LRG catalog over 
the SDSS Data Release 4 footprint ( |Collister et al.|2007 '). We 
remove sources classified as stars in the SDSS catalog, but do 
not use the star/galaxy classifier in the MegaZ catalog. Con- 
tamination by stars is estimated by Collister et al. (2007 ) to 
be 5%. We extended the catalog with the additional area from 
Data Release 6. Redshifts for the new area were estimated 
by a nearest neighbor match with the ugriz photometry. We 
estimate that this procedure smooths the redshift distribution 
by (Tz ~ .003 and has little effect on the overall redshift un- 
certainty, which is at the = .05 level. 

The CMB temperature map we used is an inverse- 
variance weighted combination of the Wilkinson Microwave 
Anisot wpy Prob e (WMAP) 5-year Q, V and W frequency 
maps (Hinshaw et al. 2008 ), with the foreground galactic- 
emission maps subtracted from each. Regions within the ex- 
tended temperature analysis mask (KQ75), which is a con- 
servative Galactic and point source mask, are left out of the 
analysis. The maps are pixelized in Healpix format (Gorski et 
al. 2005 ) at 7 arcminute resolution, which oversamples the 30- 
arcminute full-width, half-max beam. In excel lent agreement 
with previous results ( [Giannantonio et aL 2008j ), we measured 
a cross-correlation amplitude between our two data sets on 1° 
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Fig. 1. — Stacked regions on the CMB corresponding to supervoid and supercluster structures identified in the SDSS LRG catalog. We 
averaged CMB cut-outs around 50 supervoids (left) and 50 superclusters (center), and the combined sample (right). The cut-outs are rotated, 
to align each structure's major axis with the vertical direction. Our statistical analysis uses the raw images, but for this figure we smooth 
them with a Gaussian kernel with FWHM 1.4°. Hot and cold spots appear in the cluster and void stacks, respectively, with a characteristic 
radius of 4°, corresponding to spatial scales of 100 h^^Mpc. The inner circle (4° radius) and equal-area outer ring mark the extent of the 
compensated filter used in our analysis. Given the uncertainty in void and cluster orientations, small-scale features should be interpreted 
cautiously. 



scales of 0.7/^K. 

To find supervoids in the galaxy sample, we used the 
parameter-fre e, publicly available ZOBOV (ZOnes Bordering 
On Voidness; Neyrinck] |2008 ) algorithm. For each galaxy, 
ZOBOV estimates the density and set of neighbors using the 
parameter- free Voronoi tessellation ( |Okabe et a l.|2000||van de| 
[Weygaert & Sch aap 2007 ). Then, around each density mini- 
mum, ZOBOV finds density depressions, i.e. voids. We used 
VOBOZ ( [Neyrinck, Gnedin & Hamilton|2005| l to detect clus- 
ters, the same algorithm applied to the inverse of the density. 

In 2D, if density were represented as height, the density de- 
pressions ZOBOV finds would correspond to catchment basins 
(e.g. Platen, van de Weygaert & Jones 2007) . Large voids 
can include multiple depressions, joined together to form a 
most-probable extent. This requires judging the significance 
of a depression; for this, we use its density contrast, compar- 
ing against density contrasts of voids from a uniform Poisson 
point sample. Most of the voids and clusters in our catalog 
consist of single depressions. 

We estimated the density of the galaxy sample in 3D, con- 
verting red shift to distance according to WMAP5 ( [Komatsti ] 
et al.|2 008 ) cosmological parameters. To correct for the vari- 



able selection function, we normalized the galaxy densities to 
have the same mean in 100 equally spaced distance bins. This 
also removes almost all dependence on the redshift-distance 
mapping that the galaxy densities might have. We took many 
steps to ensure that survey boundaries and holes did not af- 
fect the structures we detected. We put a 1° buffer of galax- 
ies (sampled at thrice the mean density) around the survey 
footprint, and put buffer galaxies with maximum separation 
1° from each other in front of and behind the dataset. Any 
real galaxies with Voronoi neighbors within a buffer were not 
used to find structures. We handled survey holes (caused by 
blight stars, etc.) by filling them with random fake galaxies 
at the mean density. The hole galaxies comprise about 1/300 
of the galaxies used to find voids and clusters. From the final 
cluster and void lists, we discarded any structures that over- 
lapped LRG survey holes by >10%, that were < 2.5° (the 
stripe width) from the footprint boundary, that were centered 



on a WMAP point source, or that otherwise fell outside the 
boundaries of the WMAP mask. 

We found 63 1 voids and 2836 clusters above a 2a signifi- 
cance level, evaluated by comparing their density contrasts to 
those of voids and clusters in a uniform Poisson point sample. 
There are so many structures because of the high sensitivity 
of the Voronoi tessellation. Most of them are spurious, arising 
from discreteness noise. We used only the highest-density- 
contrast structures in our analysis; we discuss the size of our 
sample below. 

We defined the centers of structures by averaging the posi- 
tions of member galaxies, weighting by the Voronoi volume in 
the case of voids. The mean radius of voids, defined as the av- 
erage distance of member galaxies from the center, was 2.0°; 
for clusters, the mean radius was 0.5°. The average maximum 
distance between void galaxies and centers was 4.0°; for clus- 
ters, it was 1.1°. For each structure, an orientation and ellip- 
ticity is measured using the moments of the member galaxies, 
though it is not expected that this morphological information 
is significant, given the galaxy sparseness. 

3. IMPRINTS ON THE CMB 

Figure[T]shows a stack image built by averaging the regions 
on the CMB surrounding each object. The CMB stack cor- 
responding to supervoids shows a cold spot of -11.3/xK with 
3.7(7 significance, while that corresponding to superclusters 
shows a hot spot of 7.9/iK with 2.6(t significance, assessed 
in the same way as for the combined signal, described below. 
Figure |2] shows a histogram of the signals from each void and 
cluster 

To assess the significance of our detection, we averaged 
the negative of the supervoid image with the supercluster im- 
age, expecting that the voids would produce an opposite sig- 
nal from the clusters. We used a top-hat compensated filter 
to measure the fluctuations, averaging the mean temperature 
within 4° of the center, and then subtracting the mean tem- 
perature in a ring of the same area around it. This filter is 
insensitive to CMB fluctuations on scales larger than the ob- 
ject detected; for an uncompensated filter, these fluctuations 
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Fig. 2. — Histograms of the signals of the 50 highest-significance 
supervoids and superclusters used for our detection, measured in our 
compensated filter The vertical dotted lines are the means of each 
distribution, at -1 1.3±3.1/iK( voids) and 7.9±3.1/iK(clusters). 

would constitute a significant source of noise. 

What is the likelihood that our results are due to random 
fluctuations? To estimate that, we performed two sets of 1000 
Monte Carlo simulations. First, we generated random posi- 
tions of voids and clusters within the survey and stacked the 
corresponding areas of the actual CMB map. This models the 
errors given the observed CMB sky and foreground subtrac- 
tion, but might not properly account for any covariance due 
to the actual configuration of voids and clusters. Second, we 
generated model CMB skies smoothed to WMAP resolution 
and repeated our analysis on these with the actual void and 
cluster configurations observed in the catalogs. We find that 
these two approaches produce identical distributions consis- 
tent with Gaussians, and with standard deviations within 2% 
of each other. The hypothesis that the signal arose from ran- 
dom fluctuations is excluded at the 4.4f7 level, a 1:200,000 
chance. Our final mean signal with errors is 9.6 ± 2.2/iK. 

We note that the radii of the structures found by 
ZOBOV/VOBOZ are typically less than 4°. One possible rea- 
son is that the algorithm could be conservative in defining 
edges in the face of significant discreteness noise. The de- 
tected structures could just be the tips of larger hills and val- 
leys in the potential. The stacked signal is also likely smeared 
somewhat from noise in determining the structures' centers. 

Our procedure does have two parameters: the number of 
objects used to generate the stacked image, and the filter size 
used to assess the hot and cold spots' significance. We used 
the same number of voids and clusters for simplicity. Density- 
contrast thresholds of 4, 3 and 2a give 7, 51, and 631 voids, 
respectively. With too few structures, the measurement would 
be swamped by CMB fluctuations (with a standard deviation 
of 22/iK in our filter). With too many, structures would be 
introduced that have dubious physic ality. We used the 50 
objects with the highest density contrast (a cut at ^ Scr for 
voids, and ~ 3.3(7 for clusters) to roughly balance these ef- 
fects. Stacking 70 voids and clusters gives a signal of 2.8ct; 
with 30, the signal remains above 4a. These results appear 
in Table [T] It is almost certain that some number of objects 
would give a higher significance than 50 gives, but we did not 
search this parameter space, to simplify the interpretation. 

As mentioned earlier, the ISW signal is expected to peak at 
a radius of ^ 4° according to theory; this is also confirmed 
by the visual appearance of the images. Changing the filter 
radius between 3-5° results in various detection significances 
of approximately 4a; these results are listed in Table [T] In a 
strict Bayesian sense, even inspecting the image by eye prior 
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to statistical analysis complicates the interpretation due to a 
posteriori bias issues. This would be difficult to quantify, but 
its effect should be small because the signal was robust in the 
few cases we checked. 

There are systematic effects from foreground contamina- 
tion that, in principle, can mimic the ISW signal. Dust emis- 
sion from the Milky Way is bright at microwave frequencies 
and is correlated with the dust extinction correction used in 
the galaxy catalog. Extragalactic radio sources correlated 
with luminous red galaxies could also potentially contribute 
to a false signal at microwave frequencies. To check that we 
are not significantly contaminated by microwave sources, we 
repeated our analysis on the individual frequency maps us- 
ing the KQ75 mask as in our combined analysis, but without 
subtracting the foreground template maps. The mean ampli- 
tudes of the void signal in the Q,V and W bands were —10.6, 
— 11.1 and —11.1 /iK; the mean amplitudes of the cluster sig- 
nal were 7.8, 7.9 and 7.7 /iK; the error on each of these means 
is 3.1 /iK. These results agree with our measurement made on 
the combined map, and demonstrate that there is no signifi- 
cant frequency dependence of the signal. Moreover, the void 
signal is expected to be less sensitive to foreground contami- 
nation. 

4. DISCUSSION 

We have measured a 4a temperature deviation on the CMB 
due to supervoids and superclusters at z ~ 0.5. The most 
likely explanation for this is that we detect the ISW effect. 
The linear ISW effect vanishes in a flat universe without dark 
energy, and the higher-order ISW contribution is expected to 
be significantly lower than the ISW in ACDM ( [Seljak [1996 



|Tuluie et al.|1996[|Crittenden & Turok|1996| l. The consensus 
in the literature is that detecting the ISW effect signals the 
presence of dark energy in a flat Universe. 

To estimate the expected effect from ISW in a ACDM uni- 
verse, we measured the si gnal that the M illennium cosmolog- 
ical A^-body simulation ( jSpringel et al.|[2005) produces. We 
ray-traced through the simulation, summing up the change 
in potential that a photon would experience passing through 
the 500 h^^Mpc box in each Cartesian direction. In this vol- 
ume, which is large enough for 1 or 2 supervoids and super- 
clusters, we checked that the linear part of the ISW signal 
through the box dominates over higher-order effects. Center- 
ing a 100 h^^Mpc aperture around the maximum ISW sig- 
nal in the Millennium volume gives 4.2/iK, ^ 2a lower than 
what we observed in our CMB stack. Though we only ex- 
pect these numbers to agree to within an order of magnitude, 
we note that most previous ISW measurements are also some- 
what higher than the predicted signal in a ACDM cosmology 
( |Ho et al. 2008). While more theoretical studies are needed to 
turn our detection into precision constraints on cosmological 
parameters, we interpret our image as the ISW effect on the 
CMB caused by the decaying of potentials in an accelerating 
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universe with dark energy. 

Previous works used the two-point cross-correlation func- 
tion of 2D projected galaxy density maps with the CMB to de- 
tect the ISW effect, reaching a significance of 2 to 2.5(t for the 
galaxy sample we analyzed ( |Ho et al. 2008; Giannantonio et| 
[al.,200 8 ). Several factors likely contribute to the higher signif- 
icance of our measurement. First, we analyze only superstruc- 
tures that should be strong ISW sources. Second, we use 3D 
information to identify them. The 2D projected galaxy den- 
sity is typically not extremal at the superstructures' locations; 
thus, the cross-correlation function is not especially sensitive 
to their contributions. Third, galaxy autocorrelations directly 
contribute to the noise for the cross-correlation function, but 
not for our method. 

Our detection makes it more plausible that low-redshift 
supervoids and superclusters explain anomalies observed 
on the C MB (Rakic eraL]|2006l [Rudnick et al.|[2007l |hioue| 
|& Silk|[2007t iMaturi et al.| |2007| l. At low to moderate 
significance, these f eatures include a 5° 70/iK cold spot 
( [Vielva et al. 2004 j ), the North-South power asymmetry, 
the low quadrupole moment, and the alignment of low 
multipoles (Huterer 2006 ). Additionally, /„!, a measure of 
non-Gaussianity on the CMB, has been estimated to be pos- 
itive at low significan ce in WMAP ( |Yadav & Wandelt|2008' 



For supplementary information, i ncludi ng th e void and 
cluster catalogs, see Granett et al. ( 2008^, and http:// 



[Komatsu et al. 20081. This indicates a CMB temperature 



distribution that is slightly skewed toward low temperatures, 
as predicted by a small nonlinear ISW effect that enhances 
supervoid signals over superclusters ( Tomita & Inoue|2007| l. 
We indeed find somewhat stronger cold spots, and although 
the difference is not statistically significant, its consistency 
with the above picture is intriguing. 
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